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Abstract 



Recent proofs of the convergence of the linear delta expansion in zero and in one 
dimensions have been limited to the analogue of the vacuum generating functional in 
field theory. In zero dimensions it was shown that with an appropriate, TV-dependent, 
choice of an optimizing parameter A, which is an important feature of the method, 
the sequence of approximants Zj^ tends to Z with an error proportional to e~ cN . In 
the present paper we establish the convergence of the linear delta expansion for the 
connected vacuum function W = InZ. We show that with the same choice of A the 
corresponding sequence Wn tends to W with an error proportional to e~ cv ^ . The 
rate of convergence of the latter sequence is governed by the positions of the zeros of 
Zn- 

PACS numbers: ll.10.Jj, 12.38.Cy, 11.15.Tk. 
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1. Introduction 

In a previous paper Q it was proved that in zero- dimensional space-time the optimized 
linear delta expansion can completely cure the problems inherent in conventional pertur- 
bation theory. A conventional perturbation expansion is a formal series in powers of the 
coupling constant. Such a power series typically has a zero radius of convergence and 
sometimes is not even useful as an asymptotic series because it is not Borel summable. In 
contrast, the optimized linear delta expansion produces a sequence of approximants which 
converge rapidly to the exact answer, with an error Rn that decreases exponentially with 
the order N: R N ~ exp(— cN). The proof given in Ref. JI| was subsequently extended 
to one-dimensional space-time (the quantum anharmonic oscillator), where it was shown 
that for the finite-temperature partition function Rn ~ exp(— cN 2 ^ 3 ). The technique can 
also be extended to establish a proof of convergence for cutoff y>\ 3 theories (with either 
sign of the squared mass) in finite volume. 

However, the proofs of convergence given in Refs. IJ and were limited to the partition 
function (vacuum-vacuum amplitude) Z, which represents the sum of all vacuum graphs, 
disconnected as well as connected. In quantum field theory the crucial quantity to compute 
is W, the logarithm of the vacuum-vacuum function, which represents the sum of the 
connected vacuum graphs only. It is nontrivial to show that the optimized linear delta 
expansion converges for W because Zn, the iVtli approximation to Z, has zeros in the 
complex-<5 plane at which In Z^ is singular. The presence of such zeros could interfere with 
the convergence of the sequence Wn- It is the purpose of this paper to show that in zero 
dimensions the delta expansion for W does in fact converge, in spite of the zeros of Z^. 

The specific model we consider here is the zero-dimensional analogue of the Euclidean 
functional integral for a ip 4 quantum field theory, which in this case amounts to the one- 
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dimensional integral 

Z = I dxe- 9X "-^ x \ (1) 

J— oo 

where g is the coupling constant and \i the mass. In order to perform a weak-coupling 
expansion of Z in powers of g it is necessary that the mass parameter /z 2 be positive; a 
weak-coupling expansion does not exist otherwise. However, in the linear delta expansion 
the value of fi 2 is immaterial. Thus, for simplicity we restrict our attention to the massless 
case /i = 0; the analysis of the massive case (with either sign of fi 2 ) does not differ in any 
significant way. 

The linear delta expansion || has features in common with a number of previous ap- 



proaches P-pil to improving on the convergence of ordinary perturbation theory. It in- 
volves the introduction of an artificial parameter 5 which does not appear in the original 
problem and which interpolates linearly between the theory we hope to solve, with action 
S, and another soluble theory, with action So. The interpolating action S(5) is defined as 

S(S) = A(l - S)S + 6S, (2) 

so that S(0) = So and S(l) = S. 

Any desired quantity is evaluated as a perturbation series in powers of 5, which is 
then set equal to 1 at the end of the calculation. At 5 = 1 the theory defined by S(5) is 
independent of the value of A. However, to any finite order N in 5 there is a residual A 
dependence, and the choice of A is in fact crucial to the convergence of the delta expansion. 
Indeed, if A were taken to be a constant independent of the order N, the delta expansion 
would have a zero radius of convergence, just as in ordinary perturbation theory. However, 



it was proved in Ref. QTJ] that if A is chosen as yaN, the sequence of approximants ^(A^) 
converges to Z. The numerical value of a is given at the beginning of Section 3. The scaling 
A = cN 2 ^ 3 was shown in Ref. 0] to guarantee the convergence of the corresponding sequence 
for the finite-temperature partition function of the anharmonic oscillator. However, the 



proof did not extend to zero temperature because the limits (3 — * oo and N — > oo are not 
interchangeable. 

In the context of field theory it is natural to work with the logarithm of the partition 
function, i.e. with connected diagrams. It is therefore desirable to extend the proof of 
convergence to W = In Z. To do so it is necessary to determine the location of the zeros 
of Zn in the complex-5 plane in the limit of large N. This determination is performed 
asymptotically by a steepest-descent evaluation of the integral representing Z(8) and an 
asymptotic analysis of the behavior of the remainder as a function of complex S. 

In Section 2 we explain how the zeros of Z N affect the convergence of the 5 expansion for 
W. The asymptotic analysis of the location of these zeros is given in Section 3. Finally, in 
Section 4 we summarize our results and discuss possible extensions of the analysis presented 
here to higher dimensions. 
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2. Relevance of the Zeros of Zjy 

In this paper we are investigating the convergence of the delta expansion for W = In Z. This 
involves computing the iVth partial sum Wn of the Taylor series in 8 of W{8) = \nZ(8), 
where 

/OO 
-oo 

in which we have taken g = 1 without loss of generality. Then 8 is set equal to 1 and A 
chosen in some appropriate fashion as a function of N. 

It was shown in Ref. [|TJ that with the appropriate scaling of A the sequence of approx- 
imants Z N {8) evaluated at 8 = 1 tends to the exact result. Consequently the sequence 
In Zn also tends to W. However, this does not constitute the systematic expansion of W 
in powers of 8 that we seek. That is to say, Wn, the sum of the first N terms of the Taylor 
expansion of InZ is not the same as lnZ/v, the logarithm of the sum of the first N terms 
of the Taylor expansion of Z. 

To examine the convergence of Wn to W we will make use of a slight generalization of 
some identities introduced in Ref. [0. Given a function F(S), the iVth partial sum of its 
Taylor expansion evaluated at 8 = 1 can be represented as the contour integral 

where Cq is a closed anticlockwise contour encircling the origin but not the point z = 1. 
The quantity F(l) itself can be represented as 

where C\ is a closed contour encircling the point z = 1 but not the origin. Subtracting 
Eq. (4) from (5) gives the following general integral representation for the remainder = 

F(i)-F N (iy. 

^ 1 f dz 1 T _ j v 
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where the contour Coi encircles both z = and z = 1 in an anticlockwise direction, as 
shown in Fig. 1. 

Now let us apply these identities to the function F(5) = \nZ N {5) evaluated at 5 = 1 in 
order to obtain a bound on the remainder, which in this case we denote by 

TZ N = In Z N - (In Z N ) N . (7) 

Expressing Z^(z) in terms of its roots: 

N / z\ 

Z N (z) = Z N (0)UU--) (8) 

so that 

N / z\ 

lnZ N (z) = lnZ N (0) + - -) . (9) 

Thus, 

^-^El^rr^M— )■ do) 

Note that constant terms such as lnZ^O) do not contribute to the integral around the 
contour Cqi. 

In the derivation of Eq. (10) we are assuming that there are no singularities of the 
integrand inside the contour other than those at z = and z — 1. We will verify this 
assumption in Section 3. Let us now expand the contour Coi by pushing it outward in all 
directions. In so doing we encounter the logarithmic branch points emanating from the 
roots z r . We take the branch cuts to lie along straight lines radiating directly outwards, as 
shown in Fig. 2. Thus, the shifted contour wraps around these branch cuts and the integral 
is then given by the sum of the discontinuities across each branch cut: 

N r°° dz 1 



r= l J z r Z^^ Z 1 
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By parametrizing the integration variable z along each branch cut by pe ldr it is easy to 
establish the bound 

^ 111 

^^l^MNW (12) 

where we have used the inequality \z — 1] > \z r — 1| along each radial cut. This latter 
inequality is valid provided that \z r \ > 1, which will be established in the next section. If 
Zmin denotes the root having the smallest modulus then Eq. (12) may be replaced by the 
simpler inequality 

1^1 < r 1 —\T^Tn ■ (13) 
Since \z min \ > 1, we may conclude from Eq. (13) that the remainder H N vanishes for 
large N. The rate at which 7Zn tends to is crucially dependent on the behavior of z min 
as a function of N, which is the subject of the next section. 
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3. Saddle-Point Determination of the Smallest Zero of Zn 

In Ref. it was shown that the sequence converges to Z if the parameter A is chosen as 



A = \n = y/ctN, where the numerical value a = 1.3254 ... is obtained from a = 2/ sinh (3 
and (3 satisfies the transcendental equation (3 = coth (3. We adopt the same choice for the 
parameter A here. In this case the series 

N 

Z N (z) = J2 C nZ n (14) 
n=0 

looks rather simple: apart from the last term, which is small and negative, the coefficients 
are all positive and monotonically decreasing, looking roughly like a geometric series. Un- 
fortunately, it is not easy to determine the location of the roots of a polynomial from 
the knowledge of the coefficients c n . As an example, consider the simple- looking case 

2 

c n = e~ 7 ™ . For small values of 7 the roots lie on a circle centered at the origin in the 
complex-z plane. However, as 7 increases past a critical value two pairs of complex con- 
jugate roots break away from the circle. Additional roots eventually break away from the 
circle as 7 increases through a whole sequence of critical values. This phenomenon is de- 
scribed in Ref. |15||. To our knowledge there is no simple analytic way to determine these 



critical values or, indeed, the positions of the roots. 

Interestingly, the configuration of the roots of Z^ shares many of the characteristics of 
this simple model. For large odd N all but five of the roots lie almost exactly on a circle; 
of the remaining roots, one complex-conjugate pair lies inside the circle, another conjugate 
pair lies outside, and there is a single root located far away on the positive-real axis (see 
Fig. 3 for the case N = 27). When N is even the only qualitative difference is that there is 
no real root. As shown in Eq. (13), it is the position of the pair of roots nearest the origin 
which determines whether or not the sequence Wn converges to W. 

As remarked above, it is very difficult to find the positions of the roots directly from 
the coefficients. Thus, we have adopted the alternative strategy of expressing Z^(z) as the 
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difference 

Z N {z) = Z{z)-R N {z) (15) 

and estimating each term on the right side of (15) by asymptotic methods. Note that while 
Z(z) does not depend explicitly on the large parameter N, when z ^ 1 it does involve Ajv, 
which is a large parameter. It is through this dependence that we are able to estimate the 
integral representation for Z(z) by steepest-descent analysis. The result of this analysis is 
that 

/ 3tT \ 1 / 2 / 1 \ 

= 1+ b?) + ° i» • (16) 



From this asymptotic relation we conclude that TZn tends to zero like exp(— ^j3iiN/a), 
thus proving that Wn converges to W. 

The asymptotic analysis of the two terms Rn{z) and Z{z) on the right side of Eq. (15) 
follows in Subsections 3A and 3B, respectively. 

A. Asymptotic Bound on Rn(z) 

The starting point for our analysis is the identity for Qn(v) = e~ y {e y }N established in 
Ref. 1: 

±Q N = -y N e~y/N\ . (17) 
dy 

This identity is to be applied under the x integration with y = z(Xx 2 — x A ). That is, 



Rn(z) 



2/ dxe- Xx+y \l-e N (y)} 
Jo 



Let us first deal with the case Re z > 0, where the x integration contour can be 
maintained along the real axis. For complex z in this region we integrate Eq. (17) along 
the ray y' = yt, with < t < 1, to obtain 

8*(V) = 1- f y dy'y' N e-y'/N\ , (19) 



o 
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using 6at(0) = I. 

The x integration in Eq. (18) splits naturally into the two ranges, < x < vA and 
\/~X < x. In the first case y = \y\e ld , where 9 = arg(z). The contribution to the remainder 
from this range is therefore 

An{z) = JL j^ dx e -^y J lvl du u N e -^ e i { N + x)e (20) 

Thus, 

\A N (z)\ < dx e~ Xx2+ ^ cosd du u N e-"™ 9 . (21) 

The maximum of the u integrand occurs at u = N/ cos 9 > N, whereas the upper limit 
is less than or equal to Na\z\/4. Thus, provided that \z\ < 4/a, we may bound |Ajv(£)| by 

\A N (z)\ < 4t dx e- Xx2+ M cos9 \y\ N+1 e- ycosd 



N\ Jo 

2 \z\ N+1 [ VX dxe- Xx2 (Xx 2 -xY +1 . (22) 



iV! 1 1 Jo 

Apart from the factor of l^l^" 1 " 1 , the right side of Eq. (22) is precisely A^(l), which was 
shown in Ref. to be bounded by CN 3 ^e~ Na / 2 . Thus, 

\A N (z)\ <CN 3 /*\z\ N+1 e- Na/2 . (23) 

In the second of the two regions of the x integration, the quantity Xx 2 — x 4 is negative, 
so that y = cue^ e+n \ The corresponding contribution to the remainder is then 

B N ( Z ) = A r dx e~ Xx2+ y [ lVl du ^^^(N+Die+n) (24) 
N\ JV\ Jo 

so that 

\Bn(z)\<— dx e- Xx2 -\ y \ cosd du^e"™ 9 . (25) 

In this case the u integrand is a monotone function whose maximum occurs at the upper 
limit. Thus, 

\B N (z)\ < / dx e- Xx2 - Mcose \y\ N+1 e Mcosd 

n 



= ^-\z\ N+l r dx e- Xx \x A - Xx 2 ) N+1 . (26) 

Again, the right side is B^(l) multiplied by a factor of |z| + . In Ref. -Bjv(l) was 
shown to be bounded for large N in precisely the same way as Ajy(l). 
Thus, altogether 

\Rn(z)\ <CN 3 / 4 \z\ N+1 e- Na / 2 , (27) 



provided that 



\z\ < - . (28) 
a 



Now let us consider Re z < 0. The only difference is that since cos 6* < the roles of 
A N and B N are reversed. That is, |Ajv(z)| is bounded as in Eq. (23) independent of \z\, 
while \Bn(z)\ is bounded in the same way, provided that Eq. (28) holds. Altogether, the 
conclusion for \Rn(z)\ remains the same. 

B. Steepest-Descent Evaluation of Z(z) 

For Re z > the integral in Eq. (3) is an adequate definition of the function Z(z) in the 
complex-z plane. Later on in this subsection we will examine the analytic continuation of 
Z(z) to the left-half z plane. Because A = is large we evaluate the integral by looking 
for the saddle points of the exponent (f in the integrand: 

<p(x) = -\x 2 + z(\x 2 - x 4 ) . (28) 



The three saddle points satisfying <p'(x) = are xq = and x± = ±yj\(z — l)/(2z). At 
these saddle points (p(xo) = and <p(x±) = —X 2 (z — l) 2 /(4^). 

The method of steepest descents requires that for each value of complex z we deform the 
integration path in Eq. (3) to a stationary-phase contour in the complex-x plane connecting 
the original end points at x — ±oo. We find that there are two cases to consider. When 
\z\ < 1 (region A) the stationary-phase contour passes through the saddle point xq only 
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and not through the others. In contrast, when \z\ > 1 the stationary-phase contour passes 
through all three saddle points. The appropriate integration contours for these two cases 
are illustrated in Figs. 4 and 5, respectively. In the second case there are two subregions 
in the complex-^; plane to consider, region £>, in which exp[</?(x±)] is subdominant with 
respect to exp[</?(xo)] = 1, and region C in which the reverse is true. The boundary curve 
T between subregions B and C is given by the polar equation 

2r 

r : cos# = - , (29) 

r 2 + 1 ' V ; 

where z = re td . The above regions and the boundary curves are illustrated in Fig. 6. 

In regions A and B the dominant contribution to Z^(z) comes from the saddle 
point at the origin, whose contribution is a slowly-varying (nonexponential) function of 
N: Z(z) ~ A r ~ 1 / 4 . Referring to Eq. (15), we see that a zero of Z N (z) must arise from a 
cancellation between this contribution and Rn{z), which as we have seen is bounded by 
const x exp[iV(ln \z\ — a/2)]. Thus, no zero is possible for \z\ < e Q / 2 = 1.94.... Since 
this value is strictly greater than one, the zeros in regions A and B will not affect the 
convergence of TZ n- From the numerical plot in Fig. 3 it appears that the bound on \z\ is 
saturated, giving rise to a ring of zeros. 

In region C the dominant contribution to Z N {z) comes from the saddle points x±. These 
saddle points give a nontrivial exponential dependence to Zn(z) of the form exp[iVa:(,2 — 
l) 2 / (4z)]. Since this exponent is positive the necessary cancellation required for a zero can 
only occur at still larger values of \z\. This accounts for the conjugate pair of complex 
zeros lying outside the ring in Fig. 3 and the large real root. None of these zeros affects 
the convergence of TZ^. 

Finally, we examine the boundary curve T separating regions B and C. On this curve 
the exponent of the saddle points at x± is the same as that of x . There thus arises the 
possibility of a cancellation between their contributions, leaving a remainder which can be 
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compensated by Rn(z) at a smaller value of \z\ inside the ring of zeros in Fig. 3. Note that 
the smallest zero marked in Fig. 6 lies almost exactly on T, indicating that this is indeed the 
mechanism that produces the zero of smallest absolute value. In order for this cancellation 
to occur the two saddle-point contributions must be completely out of phase. The phase of 
each contribution is of course uniquely determined by the direction of the stationary-phase 
path going through the saddle point. As can be seen from Fig. 5, the stationary-phase path 
directions differ by tt/2. Thus, the required condition for cancellation is 

Im (p(x±) = tt/2 + (2n + 1)tt , (30) 

where n is any integer. 

This equation is to be solved in conjunction with Eq. (29). The condition \z\ > 1 
excludes negative values of n and in fact the smallest value of the modulus of the root z is 
obtained for n = 0. With n = the simultaneous solution of Eqs. (29) and (30) is 

z = a + Va 2 - 1 , (31) 

where a = 1 + 3iri/(Na). The accuracy of our asymptotic analysis improves with increasing 
N. For N = 27 this saddle-point analysis predicts the positions of the smallest roots quite 
accurately: z = 1.481 . . . ± (0.811 . . to be compared with the actual numerical values 
z = 1.443 . . . ± (0.835 . . as shown in Fig. 6. For larger values of N the roots approach 
their asymptotic values with an error that behaves like 1 / y/~N. 

From Eq. (31) we can determine the behavior for large N of the modulus of the smallest 
root z min . This is given in Eq. (16), which we have checked numerically by performing a fit 
to a series in inverse powers of y/N up to iV = 59. The crucial feature of Eq. (16) is that 
|z m in| approaches 1 from above sufficiently slowly that TZn, the difference between Wn and 
W, tends to zero. While the convergence is not as rapid as that of the sequence Z N , which 
converges like exp(— aN/2), it still converges like an exponential: exp(— ^j3nN/a). 

14 



Up till now we have not considered the case Re z < 0. For this case the integral 
representation for Z(z) is no longer valid as it stands because it is divergent. As z is 
rotated into the left-half complex plane the endpoints of the integration contour must be 
rotated in the opposite direction (and at one quarter of the rate) in order for the integral 
continue to to exist. A description of this analytic continuation procedure may be found in 



Ref. |L6[ for the case of boundary conditions on differential equation eigenvalue problems. 
Once the integral representation for Z(z) has been continued to negative values of Re z it 
may be subjected to the same saddle-point analysis as above. For this case, there are just 
two regions, \z\ < 1 in which the stationary-phase contour passes only through the saddle 
point xo and \z\ > 1 in which the contour passes through all three saddle points, but the 
contribution from xq dominates. These two cases correspond to what happens in regions 
A and B for Re z > 0. Thus, the resulting zeros complete the ring of zeros on Fig. 3 but 
do not affect our conclusions regarding the convergence of Wn- 
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4. Discussion and Conclusions 

The arguments presented above establish the convergence and bound the remainder in an 
optimized expansion of W = In Z for the non-Gaussian integral (3), provided we choose to 
vary the optimizing parameter A with iV at large N so as to optimize the convergence of 
the partials Z N (i.e. we take A ~ y/N): 

In Z - In Z N ~ e~ Na/2 

(ignoring power prefactors). With this choice of N dependence for A, we have shown that 

\nZ N - {\nZ N ) N ~ e~^ 3TN/a . 

The remainders in the partials Wn = (In Zn)n for the connected function W are thus 
asymptotically 

w - w N ~ e -^/ 3nN/a . 

Although the numerical illustrations have been for odd N, we should emphasize that 
with A chosen as V 'aN the results for even N interpolate smoothly between those for 
odd N. The distinction originally arose [Q] in the context of the principle of minimal 
sensitivity (PMS, 0) as applied to Zjv(A). There it turns out that for odd iV there is a 
single stationary point of Zn as a function of A, whereas for even iV there is a point of 
inflection but no stationary point. The present paper marks a step away from reliance on 
the PMS philosophy. Instead we are simply choosing A in such a way that the convergence 
of the sequence Wn is guaranteed. This methodology eliminates the arbitrariness that 
often occurs in applying the PMS criterion. That is, there may well be several stationary 
points (17| , and one then has to choose between them by some further criterion which itself 
needs to be justified. In the case of Wn there are indeed several stationary points, and one 
no longer has a strict inequality like Z^ < Z to help one distinguish between them. For 
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example, for N — 19 there are two maxima and one minimum in A. The first maximum 
and minimum, illustrated in Fig. 7, are reasonably close to the exact value of W, while the 
second maximum exceeds it by some 0.4%. The value of A given by A = \faN is slightly 
lower than the position of the first maximum, and gives a better estimate of W: 0.5948757 

compared with the exact value of 0.5948753 

The convergence of the optimized expansion for W is slower than that for Z with this 
choice of A(iV). It is of course possible that a more rapid convergence (with \n(W — W^) — 
—N", | < v < 1) might be obtained with a different choice of X(N). Our main object here 
has been to provide an existence proof for a convergent procedure for the connected gen- 
erating function. An understanding of the convergence at the level of connected quantities 
is crucial in higher dimensions, where the delta expansion for the full partition function 
converges at any finite spacetime volume, but at a rate which deteriorates as the volume is 
increased. As W is linear in the volume (for large volume) a convergent procedure at any 
finite volume will be uniformly convergent (for connected quantities) as the volume cutoff 
is removed. Finally, we note that the techniques used above involve only saddle-point esti- 
mates which should generalize readily to functional integrals defining the partition functions 
for quantum mechanics or field theories. 
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Figure Captions 

Fig. 1 Contours in the complex-z plane used for the representation of F N (l), F(l), and 
R N (Eqs. (4)-(6)). 

Fig. 2 Logarithmic branch cuts radiating from the zeros of Z^(z) in the representation of 
TZ N in Eq. (11). 

Fig. 3 Location of zeros of Z^(z) in the complex plane for N = 27. The circle (of radius 
2.1) is just an empirical fit to the ring of zeros. The distant real root at z — 14.124 
is not shown. 

Fig. 4 Stationary-phase contours of the function ip(x) in Eq. (28) for a typical complex 
value of z with \z\ < 1. The original integration contour lying along the real axis 
must be distorted into the contour marked AxqB, passing through the saddle point 
xq. The other saddle points play no role in the asymptotic evaluation of the integral. 

Fig. 5 Stationary-phase contours of the function (p(x) in Eq. (28) for a typical complex 
value of z with \z\ > 1. The original integration contour lying along the real axis 
must be distorted into the contour marked Ax-Bx Cx + D, passing through the three 
saddle points Xo and x±. 

Fig. 6 Right-half z plane showing the three regions A, £>, C discussed in Subsection 3B. 
The boundary between regions B and C, labeled T, is given in Eq. (29). Also shown 
are the two smallest-modulus roots of Z N (z) for iV = 27. 

Fig. 7 Wn(X) versus A for iV = 19, showing the lower maximum and the minimum. The 
dotted line is at A = = VaN. Note that the vertical scale is highly magnified, 
with a range of about 7 x 1CT 5 . 
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